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The production of the three jets in hadron-hardon collision is the first more complex process 
which allow us to define a branch of variables in order to do more precise measurement of the 
strong coupling and the parton distribution function of the proton. This process is also suitable for 
studying the geometrical properties of the hadronic final state at hadron colliders. This requires 
next-to- leading order prediction of the three-jet observables. In this paper we describe the theoretical 
formalism of such a calculation with sufficient details. We use a the dipole method to construct 
Monte Carlo program for calculating three-jet observables at next-to-leading order accuracy. We 
present a theoretical prediction for inclusive and exclusive cross section and for some relevant event 
shape variables like transverse thrust, transverse jet broadening and Eta variable. 
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The main difficulty of the next-to-leading order cal- 
culations is the presence of the singularities. In general, 
when evaluating higher order QCD cross section, one has 
to consider real-emission contributions and virtual cor- 
rections and one has to deal with different type of the sin- 
gularities. The ultraviolet singularities presented in the 
virtual contributions and they are removed by the renor- 
malization. The infrared singularities presented both in 
the virtual and real contribution but the sum of them is 
finite. In the last few years the theoretical developments 
made possible the next-to-leading order calculation for 
three-jet quantities. There are several general methods 
available for the cancellation of the infrared divergences 
that can be used for setting up a Monte Carlo program 
0, S 0| • In computing the NLO correction we use the 
dipole formalism of Catani-Seymour '9] that we modified 
slightly in order to have a better control on the numeri- 
cal calculation. This scheme is discussed in details in the 
next section. 

The advantages of using the dipole method are the fol- 
lowings: i) no approximation is made; ii) the exact phase 
space factorization allows full control over the efficient 
generation of the phase space; iii) neither the use of color 
ordered subamplitudes, nor symmetrization, nor partial 
fractioning of the matrix elements is required; iv) Lorentz 
invariance is maintained, therefore, the switch between 
various frames can be achieved by simply transforming 
the momenta; v) the use of crossing functions is avoided; 
vi) it can be implemented in an actual program in a fully 
process independent way. 

The important theoretical development that made pos- 
sible the three-jet calculation was that the relevant one- 
loop amplitudes for the relevant subprocesses became 
available. For the 99999 E3i ~^ Q'1999 subpro- 
cesses the amplitudes were calculated by Bern, Dixon 
and Kosower [IJ and for the — *■ qqQQ9 subprocess it 
was given by Kunszt, Signer and Trocsanyi |12|. The 

relevant six par ton tree-level matrix elements are also 

available El 111 IH El . 

1nagyz@physics.uoregon.edu1|http : //www ■ cpt ■ dur ■ ac . uk/~nagyz/nio++| In this paper we give Sufficient details of Our work and 



I. INTRODUCTION 

In high energy hadron-hadron collision the processes 
with purely hadronic final state is one of the clearest 
processes to test the Quantum Chromodynamics (QCD) 
an measure its parameters simultaneously. According to 
the perturbative QCD the cross section of this processes 
is a convolution of a long- and a short distance part. 

The long distance dependent part is the parton distri- 
bution function of the incoming hadrons. This function 
is universal, process independent and we can measure it 
in any basics process for example deeply inelastic scat- 
tering. 

The short distance part is the partonic cross section 
that can be calculated in perturbative QCD as a function 
of a single parameter, the strong coupling (as)- The main 
advantage of the hadronic final states at hadron-hadron 
colliders is that it is possible to measure the strong cou- 
pling and the PDF function simultaneously. 

The production of three jets is the first process which 
can provide complex final state to define a branch of jet 
observables in order to be able to do more precise mea- 
surement of the strong coupling and give better determi- 
nation of the parton distribution function. On the other 
hand this process allow us to do more advanced studies 
of the hadronic final states by measuring its geometri- 
cal properties. In order to be able to make quantitative 
prediction it is essential to perform the computations at 
least next-to-leading (NLO) accuracy. In hadron collision 
the most easily calculated one- and two-jet cross sections 
have so far been calculated at NLO level EiI2]. At next- 
to-leading level some inclusive threedet observables were 
calculated by Giele and Kilgore H 0| and by the author 
1^. Furthermore, Z. Trocsanyi also calculated the three- 
jet cross section in effective parton distribution function 
approximation 
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present several new results for next-to-leading order pre- 
diction of three-jet observables in hadron-hadron colli- 
sion that were not published before and they could be 
interesting in an experimental analysis. In the Sec. ^ 
we give details of the analytic and numeric calculation 
and describe the structure of our result. In the Sec. IfffI 
we present complete NLO predictions for three-jet inclu- 
sive cross sections and their energy fraction distribution 
(Dalitz variables) using the inclusive fcj_-algorithm [T^ 
and the midcone algorithm j^^. We define the thrust 
and jet broadening event shapes variables in the trans- 
verse plane. We show some subjet result using the exclu- 
sive /c^-algorithm "IJl. Sec. IIVI contains our conclusion. 



II. DETAILS OF THE CALCULATION 



A. The method 

The next-to-leading order cross section for the process 
with two initial state hadrons is the convolution of the 
parton density functions of the incoming hadrons and the 
hard scattering cross section, 

ct(p,p)=X! / dvfa/HiV.fJ'F) / dfj f^/HiV.fJ'F) 
,. I. JO Jo 



where f^/j^ is the density of the partons of type i in the 
incoming hadron H at the factorization scale. The 
corresponding parton level cross sections are 

a^,"{p,p) ^ [ da^,{p,p) = / dr(3)(|A^(3)|2)Ff , (2) 



and the next-to-leading order correction is sum of three 
integrals 
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{P,P) 



dOab{P,P) 



dcrIbiP:P)+ I da'^^{p,p) , (3) 



where da^ and da^ are the real and virtual contributions 
to the cross section. The contribution da'^ represents the 
collinear subtraction counter term. The pole structure of 
this term is well defined but the finite part is factor- 
ization scale dependent. The actual form of the finite 
part in the expression of the da'~^ defines the factoriza- 
tion scheme. In our calculations we chose the commonly 
used MS subtraction scheme. The parton density func- 
tions are also scheme dependent, so that this dependence 
cancels in the hadronic cross section of Eq. 

The three integrals on the right hand side of Eq. Q 
are separately divergent but their sum is finite provided 
by the jet function Fj"''' defines an infrared safe quantity, 
which formally means that 
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in any case where the four-parton and three-parton con- 
figurations are kinematically degenerate. The presence of 
the singularities means that the separate pieces have to 
be regularized and the divergences have to be cancelled. 
We use the dimensional regularization in d = 4 — 2e di- 
mensions in which case the divergences are replaced dou- 
ble poles 1/e^ and single poles 1/e. We assume that ultra- 
violet renormalization of all Green functions to one-loop 
order has been carried out, so the divergences are infrared 
origin. In order to get the finite sum a slightly modified 
version of the Catani-Seymour f9| dipole method is used 
in our calculations. 

The reason for modifying the original dipole formalism 
is numerical. The essence of the dipole method is to de- 
fine a single subtraction term da"^ , the dipole subtraction 
term, that regularizes the divergences in the real term 
which comes form the soft and collinear regions. Thus, 
the three singular integrals in Eq. ||2J) are substituted by 
three finite ones: 

a:riP,p) = '^^bf?}iP,P) + -lUiP'P) 



dx 



cT^bmi^,xp,p)+a^,\^yix,p,xp) , (5) 



where the four-parton integral is given by 



NLO I 



x[af'°(7yp,77p)+af,^0(,7P,^p)] , (1) "-H3}^ 



iP,P)= / [d<7ab{P,P)e=0 - da^f,{p,p),=o] 



(6) 



We have two three-parton contribution to the NLO cor- 
rection. The second term on the right hand side of Eq. 

is the sum of the one-loop contribution and a Born 
term convoluted by an universal singular factor I, 



[daX,{p,p)+da^,ip,p)^l\^^^ . (7) 



The factor I contains all the e poles which come from 
the da^ and da'-' that are necessary to cancel the (equal 
and with opposite sign) poles in da^ . The last term 
in Eq. ||SJ) is a finite remainder that is left after fac- 
torization of initial-state collinear singularities into the 
non-perturbative distribution functions (parton density 
function). 



/ dx a^i^f^-. {x, xp, P) = V / dx 
JO „, Jo 



daj^,,{xp,p)<E>{P{x)+K{x)r 



, (8) 



(4) 



where the x-dependent functions P and K are similar 
(but finite for e ^ 0) to the factor I. These functions 
arc universal, that is, they are independent of scattering 
process and of the jet observables. 

There are several ways to define the da^ dipole sub- 
traction term, but all must lead to the same finite next- 
to-leading order correction. In this program a slight mod- 
ification of the Catani-Seymour subtraction term was im- 
plemented by defining the dipole term as a function of a 
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parameter a £ (0, 1] which controls the volume of the don't want to give the precise definition of every vari- 



dipole phase space. In the e^e~ annihilation case this 
modification of the dipole formalism was discussed in 
Ref. |22| what we generalize for the hadron-hadron colli- 
sion. 

In this subsection we recall only those dipole factoriza- 
tion formulaes that are relevant in our calculation. We 



J 



ables, functions, factors we just use the same notation of 
the original paper of the dipole method Q . 

The da^ local counter term is provided by the dipole 
factorization of the tree level squared matrix element. 
Thus we can define 



dai 



XS X! X! '^'■j,k{Pa,Pb,Pl, ■■■,Pi)F^f\pa,Pb,Pl, ■■,P^J,Pk, ■■,Pi)Q{yi]M < a) 

K, pairs ky^ij 

+ X! [^ijiPaiPbiPl, ■■■,Pi)Ff\pa,Pb,Pl,-;Pij, ..,P4)6(1 - Xy,a < Of) -f (a 5) 

+ \^'ki.Pa,Pb,Pi,--,Pi)Ff\pa,Pb,Pi,-,Pk, ..,P4)6(ui < a) + (a ^ 6) 
+ J2[T^''"-\Pa,Pb,Pi,:;P4)F^'^\pa,Pb,Pi,-,Pi)Q{n, < a) + (a ^ 6)] | 



(9) 



where dF^"') is the four-parton phase space including all 
the factors that are QCD independent, J2{4} denotes the 
sum over all the configurations with 4 partons and 5{4} 
is the Bose symmetry factor of the identical partons in 

Vf} and V^''' are the 



the final state. The V, 



dipole terms. The function F) is the jet function which 
acts over the three-parton dipole phase space. The j/y.fc, 
Xij^a, Ui, Vi are the dipole variables defined by the dipole 
factorization of the phase space. The dipole factoriza- 
tion of the phase space is exact phase space factorization 
which means there is no approximation used in the kine- 
matically degenerated regions (soft, coUinear and soft- 
collinear). We found it makes the Monte Carlo integral 
more stable because the real contribution and the sub- 
traction terms are defined in same phase space point. 

In Eq. lO the a = 1 case means the full dipole sub- 
traction and it gives back the original dipole subtrac- 
tion terms. In a computer program the large number of 
dipoles terms and complicated analytic structure of the 
expressions makes the evaluation of the subtraction terms 
rather time consuming. Using this cut dipole phase space 
we can speed up the program. The parameter a is also 
useful to check our program by varying the value of a 
and checking whether the full correction is independent 
of this parameter. 

The most serious numerical defect of the subtraction 
schemes is the missed binning. This happens when a huge 
positive weight from the real part and the corresponding 
huge negative weight form the subtraction term are filled 



r 



into different histogram bins. It is obvious this cut can 
increase the numerical stability of the Monte Carlo pro- 
gram by decreasing the size of the dipole phase space 
which reduces the chance of the missed binning. 

The introduced a phase space cut parameter requires 
to redefine the I and K(a;) fiavour kernels. The singular 
factor I is given by 



...,pm;a;e) 



1 



where the indices /, J runs over both final and initial 
state parton and the pi ■ pj dot-products are always pos- 
itive and 



V,(o,e) 



1 

+ 7- 
e 



+ 7,:+i^»(")+0(e), 



where 



K, - T- In^ a -I- 7i (a - 1 - Ina) 



(11) 



(12) 



The Ki and % constants are defined in the Eq. (5.90) 
and Eq. (7.28) in Ref. 0. The K{x,a) flavour kernel is 
given by the foUowings: 
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l-x 



+ a (5(1 — a:) 



(13) 



where Kp'^g (x) is defined by the factorization scheme. In the case of the MS scheme this functions are identically 
zero. The flavour functions K'^''{x,a) are 



K 



■"'{x^a) = P'^^[^)+p-^(x)\n^^^^ + 5^^Tl5{l~x)\n^a 



gab 



"W-X 



In 



1 - X 



-5{l-x)[ja + Kaia)--7r'Tl 



In — ^ '- - In Q(x <l-a 



"1 — xV 1 + a — X 1 — X 



(14) 



The (x, a) flavour kernels are defined by 



K-\x, a) = P^Ux) In < P,,b(x, e - 0) > lna(a;) - S^^TI 



' .m. ' 



1 — X 1 — X 



+ -6{l-x) 



bT-r\2 



In —Lr — t - ln(2 - x)Q{x > I - a) 



(15) 



where function a{x) is given by the following 



a{x) — min < 1, 



l-x 



(16) 



The definition of the Altarelli-Parisi probabilities 
{P^^g{x), Pab{x)) and splitting functions (< Pa,b{x,e) >) 
are given in Appendix C in Ref. The "/3 = 1 — a" 
prescription is defined by 



~ J fj Jp i — a; 



(17) 



B. Structure of the results 

Once the phase space integrations are carried out, we 
write the NLO jet cross section in the following form: 

X [pa,Pb,as{^l\),^J?Ii/Q'\Is^^^l'/Q'Hs\ Ci8) 

where /i/n (?7i Mj^) represents the patron distribution 
function of the incoming hadron defined at the factor- 
ization scale /^F = xpQnSi Va,b is the fraction of the 
proton momentum carried by the scattered partons Pa,b, 



Qhs is the hard scale that characterizes the parton scat- 
tering which could be Et of the jet, jet mass of the event, 
etc and i^u = xrQhs is the renormalization scale. 

In the presented leading- and next-to-leading order 
results we use the C/C-I-+ implementation |^ of the 
LHAPDF library 22] with CTEQ6 "2^ parton distri- 
bution function and with the corresponding as expres- 
sion which is included in this library. The CTEQ6 
set was fitted using the two-loop running coupling with 
as{Mzo) = 0.118. 

To ensure the correctness of the result several check 
was performed: (i) the Born level two-, three- and for-jet 
cross sections were compared to the prediction of NJETS 
[l^ program and perfect agreement was found; (ii) sin- 
gular behaviour of the real and dipole subtraction terms 
were checked numerically in randomly chosen phase space 
points; (iii) the total NLO corrections was independent 
of the a parameters that controls the dipole phase space; 
(iv) the program is based on the NLOJET+-I- program 
library |^1| that was used other already well tested pro- 
cesses. 

In the NLO calculation the difference of the real con- 
tributions and the dipole subtraction terms is still sin- 
gular but these are integrable square root singularities. 
Integrating these singularities by simple Monte Carlo in- 
tegration technique (choosing random values of the inte- 
gration variables uniformly) is not efficient way because 
the variance of the estimate of the integral is infinite and 
we are not able to estimate the statistical error of the 
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integral. To improve the stability and the convergence of 
the Monte Carlo integral the phase space is generated by 
multi-channel weighted phase space generator [23| • 



III. RESULTS 

In this section we study those jet observables which 
could be interesting in the jet studies at hadron-hadron 
colliders. We show some examples of inclusive three-jet 
cross sections, energy fraction distributions of the jets, 
event shapes variables and sub-jet rates. To present in- 
clusive jet cross section and some related jet observables 
like energy fraction distribution we use the EUis-Soper 
inclusive k±_ algorithm and the midpoint cone algo- 
rithm 01 . To do sub-jet studies and calculate some jet 
clustering algorithm related event shape variable like flip- 
ping variable we use the exclusive k± al gori thm which 
was defined by Catani et. al. in Ref. [T^. We also 
present some jet clustering independent event shape dis- 
tribution which are defined in the transverse plane (trans- 
verse thrust and transverse jet broadening). 



All jet contributes to this sum which passed the cuts. 

Calculating the transverse energy distribution of the 
leading jet and comparing the leading and the next-to- 
leading order fix order predictions we can see in the fcj^ 
case if the scales are set to xji^p — 1 the NLO cor- 
rection is about 25 — 30% as shown in Figure ^ The 
K — a^^'-' /a^'-' factor strongly depends on the scale 
choice but the NLO prediction is stable. Varying the 
scales between 0.5 and 2 the scale uncertainty of the NLO 
prediction is about 8 — 10% which is significantly smaller 
than the leading order uncertainty which is about 100% 
or more. 

Using the midcone algorithm we got a very similar re- 
sult. In this calculation the cone size is i? = 0.7 and the 
merge/split threshold parameter is / = 0.5. In this case 
the NLO correction is about 20 — 30% with the xr^f = 1 
scales. 

In Table^the NLO cross sections are tabulated. Com- 
paring the values we can see that fc_L algorithm gives 
higher cross section in the low Et region while in the 
high Et region the cone result is higher. The difference 
between the two results is less than 10% almost every- 
where. 



A. Three-jet inclusive cross sections 



B. Topology of the three-jet events 



The most important jet quantities are the inclusive jet 
cross section and its transverse distribution. In hadron- 
hadron collision for jet studies a variant of the cone algo- 
rithms was used. Unfortunately those algorithms which 
were used at the first run of the Tevatron are infrared 
and/or coUinear unsafe or only safe up to a given order in 
the perturbative calculation. It is essential that the used 
jet algorithm to be infrared safe all order otherwise the 
perturbative calculation and the perturbative prediction 
are unstable. Using an almost infrared unsafe jet algo- 
rithm we are not able to estimate the effect of the higher 
order contributions due to the variation of the renormal- 
ization and factorization scales. The other problem with 
this type of the jet algorithms is that they are logarith- 
mically sensitive to the detector resolution (energy and 
angular resolution). That is bad because in the fix or- 
der perturbative calculation it is impossible or hard to 
consider or simulate precisely these effects . 

In this subsection we use the EUis-Soper inclusive 
fc_L algorithm and the midcone iterative cone algorithm 
which are infrared and coUinear safe to all orders. In 
both cases we consider only those jets as hard final state 
jets which have at least 20GeV transverse energy and 
are in the |?7| < 2 pseudo-rapidity interval. The total jet 
transverse energy must be larger than 80GeV. The renor- 
malization and factorization scales are characterised by 
the hard scale which is the third of the total jet transverse 
momentum 



Qhs 



1 



^^t{j) 



3 ^ 



(19) 



One can study the structure of the three-jet events. 
We can define variables which test certain properties of 
the jets. These quantities could be any azimuthal corre- 
lation between the jets or the energy fraction of the three 
leading jets. 

The energy fraction of the three leading jets (leading 
in transverse energy) was measured by the CDF collabo- 
ration |2^ |23| ■ This quantity is defined in centre of mass 
frame of the three jets. The jets are ordered and labeled 
by energy in the rest frame [Ei > E2 > E^). The en- 
ergy fraction variable (usually called Dalitz variables) are 
given by 



X,. = 



2E,, 



E. 



i = 1,2,3 



(20) 



The variable Xi varies between 2/3 and 1, X2 between 
1/2 and 1, and X3 between and 2/3. There are only 
two independent variables since Xi + X2 + X^ = 2. We 
use only the first two variables. 

In the experimental analysis the iterative cone algo- 
rithm was used with an additional cut on the jets. It 
was required that the jets were well separated. This ad- 
ditional cut ensures that the cross section is calculable 
order by order in perturbation theory but the jet algo- 
rithm is still not fully infrared safe, it is almost infrared 
unsafe _25j. 

To get stable theoretical predictions in our calculations 
we use the inclusive k± and midcone jet algorithms to re- 
solve jets with the same cuts what was used in Ref. 27,5. 
All jets are required to have at least 20GeV transverse 
energy and they must lay in the jr;! < 2 pseudo-rapidity 
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FIG. 1: The fix order QCD predictions for the inclusive three-jet differential cross sections of the transverse momentum of the 
leading jet obtained using the inclusive k± and midcone algorithms. The bands indicate the theoretical uncertainty due to the 
variation of the scales xr^f between 0.5 and 2. The grey band is the leading order and the dark grey band is the next-to-leading 
order. The inset figures show the K factor and its scale dependence. The solid line represents the xr^f = 1 scale choice. 
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FIG. 2: The scale dependence of the total 3-jet cross sections. 



window. We calculate the normalised differential dis- 
tributions of the Xi and X2 variables and the double 
differential distribution of these variables. In these cal- 
culations the renormalization and factorization scales are 
set to the third of the total jet transverse energy Eq. (|19|l . 

The distributions are normalised by the total 3-jet 
cross section plotted in Figure El We can see the NLO 
correction reduced the scale dependence of perturbative 
prediction. The scale uncertainty of the NLO result is 
about 10 — 12% varying the scales on the 0.6 < xr^p < 2 
range while the LO uncertainty is much higher 38 — 40%. 
The corrections with the xji f = 1 scale choice are about 
30%. 



The double differential distributions of the variables 
Xi and X2 are plotted in Figures |21 The phase space 
would populate the available Xi, X2 region uniformly. 
Deviation from the uniform distribution shows the effect 
of the QCD dynamics. 

Taking the double differential distribution and project- 
ing on either axis the distributions of variables Xi and 
X2 can be obtained. In FigureQlthe differential distribu- 
tions of the energy fractions are plotted in the inclusive 
k±_ and midcone algorithm cases. 

Comparing the leading order result to the next-to- 
leading order result we can see that the NLO distribution 
of the leading jet is less peaked than the LO but in the 
case of the X2 distribution the NLO result is more peaked 
around 0.65. The NLO correction of the distributions is 
small. 

Comparing the two algorithms to each other we can 
say that the distributions are very similar and there is 
not too much difference. The Xi distribution of the in- 
clusive k±_ algorithm is a little bit more peaked and the 
X2 distribution is a little bit flatter than in the midcone 
case. The corresponding peaks are in the same positions. 



C. Three-jet event shapes 

We can distinguish two type of the event shapes. One 
can calculate and measure the event shapes which are 
associated to a jet algorithm. Generally we can say that 
this type of variables gives information about any geo- 
metrical property of the jets in a n-jet event. Or we can 
define event shape variables which measure any geomet- 
rical property of the event e.g thrust or jet bordering. 

Using the exclusive k±_ algorithm there is a natural 
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FIG. 4: The energy fraction distribution of the leading (Xi) 
and second leading {X2) jets. The upper figure is result with 
the inclusive k± algorithm and the lower figures shows the 
midcone result. 



way to introduce event shape variables. The prescription 
of this algorithm introduces a stopping parameter dcut 
(resolution variable), that defines the hard scale of the 




process and separates the event into a hard scattering 
part and a low-p^ part (beam- jet). 

Defining the dcut resolution variable event-by-event we 
can calculate the differential distribution of the Etn event 
shape variable. Denoting by dm the smallest resolution 
variable when the event has m hard final state jets the 
Etn event shape variable is defined by 



Etn = niax{d,„} 

m>n 



(21) 



We plotted the differential distribution of the Et3 event 
shape variable in FigureEl In this calculation the \ri\ < 2 
pseudo-rapidity cut was applied for every hard final state 
jet. The renormalization and factorization scales were 
given by ^j?,f — ^h^pQhs where the hard scale is the 
average transverse momentum of the three jets 



1 

Qhs = 



(22) 



Comparing the next-to-leading order result to the lead- 
ing order result we see that the LO prediction strongly 
depends on the scales. Varying the scale parameters xr^f 
between 0.5 and 2 the scale uncertainty of the NLO re- 
sult is about 15% in contrast to the LO scale uncertainty 
which is about 80%. Setting the scales to xr.f = 1 we 
see that the NLO correction is about 20 — 35%. 

One can define event shapes on the transverse plane. 
Important example is the transverse thrust which is de- 
fined by 



T± = max 



(23) 
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TABLE I: Next-to-leading order results for the differential 
distribution the three jet inclusive cross sections as function 
of the transverse energy of leading jet. The scales are xr^f = 
1. The errors are the statistical error of the Monte Carlo 
integration. 
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where p^^^i is the transverse component of the parton 
(hadron) momentum. The unit vector that n maximises 
the ratios of the sums is usually called thrust axis. In 
these sum only those particles are counted which are fulfil 
all selection criteria Cm- This selection must be infrared 
safe. In this calculation we required that the pseudo- 
rapidity of those particles which contribute to Tl must 
be in the [—1.1, 1.1] pseudo-rapidity window thus we have 



C 



N 



W\ < !•! 



1, 



,7V} , (24) 



where TV is the number of particles in the event. This def- 
inition of the thrust for hadron colliders gives an infrared 
safe longitudinally boost invariant quantity. 

The jet broadening is an associated quantity to the 
thrust. Using the thrust axis we can define B±^ on the 
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FIG. 5: Distribution of the Eti event shape variable at LO 
an NLO level. The bands indicate the theoretical uncertainty 
due to the variation of the scales xr^f between 0.5 and 2. The 
grey band is the leading order and the black band is the next- 
to-leading order. The inset figure shows the K factor and its 
scale dependence. The solid line represents the xr^f ~ 1 scale 
choice. 



transverse plane by the following formulae 



(25) 



In Figure El the differential distributions of the trans- 
verse thrust and the transverse jet broadening have been 
plotted at LO and NLO level. In this calculation the 
total transverse energy is larger than lOOGeV 



Ht=Y^ ET,^ > lOOGeV 



(26) 



The perturbative results logarithmically divergent at the 
edge of the phase space. This divergence occurs at T±_ ~ 
1 and T = 2/3 in the thrust distribution and at B± — 
and _B = 1/3 in the jet broadening distribution. Setting 
the scale to the total transverse energy the NLO order 
correction in the middle of the distribution where the 
effect of the logarithm is small is 30 — 35% for the thrust 
and 20 — 45% for jet broadening as it is shown in the inset 
figure. It is necessary to do all order resummation of the 
leading and next-to-leading logarithms to be able to do 
quantitative comparison of the theory to the data. 

One can define the average value of the event shape 
variables as function of the total transverse energy. This 
quantity is given bay the following ratio: 



/ da 



\dHT 



"^^^O'lTVlJr ' (27) 



where O3 could be either 1 — Tl or B± . 

Figure [7\ shows the average value of the transverse 
thrust. With the ^r.f = Ht scale choice it was found 
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FIG. 6: Distribution of 1 — T± and B_i_ event shape variable 
at LO {dashed line) and NLO {solid line) level. The scales 
were set to fJ,R,F = Ht- The inset figure shows the /("-factors. 
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FIG. 7: Average value of the 1 — T± event shape variable at 
(5 = at LO {dashed line) and NLO {solid line) level. The 
scales were set to ij.r,f = Ht- The inset figure shows the 
A'-factor. 



that the NLO correction is smah about 15%. The NLO 
correction doesn't change the shape of the distribution 
as it is shown by the K-factor plot. This result indicate 
that the NLO fix order prediction might describe the data 
itself with good accuracy. 

Figure IHl shows the average value of the transverse jet 
broadening. The CDF collaboration [2^ measured this 
distribution with i5 = choice and the data was compared 
to the leading order QCD prediction. Comparing the 
data to the NLO prediction the difference between the 
data and the theory is huge. This can happen because the 
contribution of the large logarithms from the small B± 
region is very large and in the NLO case it causes a huge 
negative effect. It is clear that the fix order calculation 
is unable to describe the data itself. The resummation of 
the large logarithms is important p9l |. 

As shown in Fig. El the large negative contributions 
come form the small B± region {B± < 0.02). Resum- 
ing the leading and the next-to-leading logarithms the 
B±da/dB± distribution is positive definite and we can 
expect that it is reasonable small in the small B± region. 
On the other hand Fig.|H|indicates that the NLO fix order 
calculation is stable prediction in the large {B± > 0.2) 
region. Assuming that the all order contribution of the 
B± < 0.02 region to the average value is small and chang- 
ing the S parameter we can make a better comparison of 
the data and the NLO prediction. In Fig. |Slwe can see 
that the {B±)q,o and the (-B_l)o.ooi NLO results suffer 
on the effect of the large logarithm but changing the 6 
parameter between 0.005 and 0.02 the NLO prediction 
is stable and hardly depends on the S parameter and 
the agreement with the data is much better. The inset 
figure shows the Ks {Bi_)^^° / {B±_)'^° factor. For 
0.005 < 5 < 0.02 the NLO correction is positive, nearly 
constant and it is maximum 30%. 
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FIG. 8: Average value of the B± event shape variable at S — 
at LO and NLO level for 5 = 0,0.001,0.005,0.01,0.02. The 
scales were set to jJ.R,F ~ Ht- The markers {diamonds) and 
the error bars represent the CDF data 28j and its statistical 
error. The inset figure shows the if-factors. 



Of course this is not a precise analysis and without 
the resummation of the leading and next-to-leading log- 
arithms the quantitative comparison of the current data 
and the theory is impossible but it can help to under- 
stand the huge discrepancy between the data and the 
NLO prediction. It also indicate that the NLO predic- 
tion can describe the data itself very well ii S > 0.005. 
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D. Sub-jet fractions 

In this subsection we study the sub-jet multipHcities of 
the events. We use the exclusive fc_L algorithm The 
jets are defined in a two-step clustering procedure. The 
first step of the algorithm identifies the low-pr scatter- 
ing fragments and includes them in the beam jets, thus 
factorizing the hard scattering subprocess. The scale of 
the hard final state jets is Ecut the stopping parameter 
of this clustering step. Defining the sub-jet resolution 
variable ycut — Qo/Ecut the second step of the algorithm 
resolves the sub-jet structure of the hard final state jets. 

The n sub-jet rate is defined by the ratio of the n sub- 
jet cross section and the total cross section as 



Rn{E 

cut ; Ucut ) 



<^n{Ecuti Vcut) 
'^tot{Ecut) 



(28) 



where the a is the total cross section defined by the sum 
of the n-jet cross section at the Ecut scale 



The fix order perturbative predictions depend on the 
unphysical renormalization and factorization scales. In 
these calculations the scales are given by event-by- event 
basis. We use the invariant mass of the two hard final 
state jets found at Et2 scale as renormalization and fac- 
torization scale 



^J■R,F = {Pi +P2) 



(35) 



FigO shows the sub-jet ratios and the sub-jet multi- 
plicity. The result is given as function of Ecut param- 
eter and the logarithm based ten of the sub-jet reso- 
lution variable. The distributions are plotted on the 
10 < Ecut < 110 GeV and -3.9 < logwiycut) < region. 
We can see the sub-jet ratios are strongly dominated by 
the large logarithms in the low Ecut and ycut regions. To 
be able to do any quantitative comparison to the experi- 
mental data in this region it is necessary to sum at least 
the leading- and next-to-leading logarithms at all order 
and match it to the fix order calculation. 



CTtot{E c 



n=2 



(7n{Ecut) 



t2 



da 



(29) 



The £"42 variable is the smallest value of the Ecut reso- 
lution variable where the event has two jets. Using the 
similar event shape variables 2/„ and y-n+i of the second 
step we can define the exclusive sub-jet cross section at 
a given Ecut scale by 



0'n{Ecut, Vcut) 



dEr 



Ecu 



da 



dEt2dy„ 
da 



(30) 



In the fix order perturbation calculation we can calcu- 
late the 3-sub-jet ratio at NLO level {R^^^ (Ecut, ycut)) 
and the 4-sub-jet ratio at LO level {R^'^ [Ecut ^ Vcut))- 
From the definition of the n-sub-jet ratios we have the 
normalisation condition 



i=2 



Rn{Ecut, ycut) 



(31) 



Using this relation the 2-sub-jet ration can be obtained 
at next-to-next-to-leading order level (NNLO) 



^2 \J^cut,ycut 



1 



^^^'^ {Ecut, Vcut) 

3^ 



^3 

R'i''{Ecut,Vcut) + 0{as-') . (32) 



Analogously to the e^e -annihilation we can define 
the sub-jet multiplicity 



00 

N{Ecut,Vcut) = '^nRn{E cut, ycut) 
n=2 



(33) 



and at up to second order 
N {Ecut, ycut) = 2 



R 



NLO 



{Ecut , ycut) 



+2R^''{Ecut,ycut) + 0{as^) . (34) 



IV. CONCLUSION 

This paper dealt with the next-to-leading order calcu- 
lation of three-jet observables in hadron-hadron collision. 
We gave a modification of the original dipole method that 
made possible to construct a Monte Carlo program for 
this process. We find this modified dipole subtraction 
method was very useful from the point of the numeri- 
cal calculation. The introduced phase space cut in the 
dipole term could reduce the evaluation time of subtrac- 
tion terms and increase the numerical stability of the 
Monte Carlo integral by reducing the probability of the 
missed binning. 

We calculated the transverse momentum distribution 
of the leading jet of the three-jet inclusive cross section 
using the inclusive k± an midcone algorithms. We found 
that the NLO correction can stabilize the theoretical pre- 
diction by reducing the renormalization and factoriza- 
tion scale dependences significantly. We also calculated 
the energy fraction distribution of the three-jet inclusive 
cross section. Comparing the leading order result to the 
next-to-leading order result we found that the NLO cor- 
rections are negligible in these analysises. 

We defined event shape variables in the transverse 
plane. We calculated the differential distribution of the 
transverse thrust and transverse jet broadening and the 
transverse energy dependence of the average values. We 
found that the size of the NLO correction is not too large 
and acceptable but at the edge of the available phase 
space the result suffers on the contribution of large log- 
arithm. In this region a next-to-leading logarithmic ap- 
proximation (NLL) matched with the fix order predic- 
tion can give reasonable result. The average value of 
the transverse jet broadening was compared to the CDF 
data and we found that the comparison of the current 
data and the theory without NLL resummation is im- 
possible. Using a rough approximation in the small B± 
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region we were able to do some qualitative comparison 
between data and the NLO QCD prediction and we found 
reasonable good agreement. 

We also calculated the subjet rates and the subjet mul- 
tiplicity using the exclusive k± algorithm. These quanti- 
ties were defined in that way to mimic e^e^ annihilation 
like situation at hadron colliders. The result is very sim- 
ilar to what we have in e"'"e~ annihilation. The subjet 
ratios are strongly dominated by the final-state large log- 
arithm. Upgrading the fix order calculation with next- 
to-leading logarithmic approximation this quantity may 



provide a sensitive measurement of the strong coupling 
at hadron colliders. 



Recently two general methods have been developed 
by Banfi, Salam and Zanderighi 29j and by Bonciani, 
Catani, Mangano and Nason |30| to carry out the next- 
to-leading logarithmic resummation for any observable. 
These methods can help us to improve our fix order cal- 
culation and to be able to make more precise comparison 
between the data and the theory. 
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